function x = sine(beta, x0, n)
%   输入参数：
%       beta - Sine 映射参数 (范围一般在 [0,1])
%       x0   - 初始值 (通常取在 (0,1))
%       n    - 迭代次数
%
%   输出参数：
%       x    - 存放每步迭代结果的向量

    x = zeros(1, n);
    x(1) = x0;
    for i = 1:n-1
        x(i+1) = beta * sin(pi * x(i));
    end
end
